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Time-dependent mode structure for Lyapunov vectors as a collective movement in 

quasi-one-dimensional systems 
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Time dependent mode structure for the Lyapunov vectors associated with the stepwise structure 
of the Lyapunov spectra and its relation to the momentum auto-correlation function are discussed 
in quasi-one-dimensional many-hard-disk systems. We demonstrate mode structures (Lyapunov 
modes) for all components of the Lyapunov vectors, which include the longitudinal and trans- 
verse components of their spatial and momentum parts, and their phase relations are specified. 
These mode structures are suggested from the form of the Lyapunov vectors corresponding to the 
zero-Lyapunov exponents. Spatial node structures of these modes are explained by the reflection 
properties of the hard-walls used in the models. Our main interest is the time-oscillating behavior 
of Lyapunov modes. It is shown that the largest time-oscillating period of the Lyapunov modes is 
twice as long as the time-oscillating period of the longitudinal momentum auto-correlation function. 
This relation is satisfied irrespective of the particle number and boundary conditions. A simple 
explanation for this relation is given based on the form of the Lyapunov vector. 

PACS numbers: 05.45.Jn, 05.45.Pq, 02.70.Ns, 05.20.Jj 



I. INTRODUCTION 

Statistical mechanics based on dynamical instability 
has drawn considerable attention in recent years. The 
dynamical instability is described as a rapid separation 
of two nearby trajectories, the so called Lyapunov vector, 
and causes the loss of memory or unpredictability of the 
dynamical system. The exponential rate of expansion or 
contraction of the magnitude of the Lyapunov vector is 
called the Lyapunov exponent, and its positivity, mean- 
ing that the system is chaotic, is a well known indica- 
tor of the dynamical instability. Many efforts have been 
devoted to connect the dynamical instability with statis- 
tical properties, like transport coefficients 0,0- Some 
works concentrated on specific effects of the dynamical 
instability in many-particle systems. Information on the 
dynamical instability in many-particle systems is given 
by the set of Lyapunov exponents (the Lyapunov spec- 
trum) and their associated Lyapunov vectors, and their 
structures have been of much interest to the study of 
chaotic properties of many-particle systems. The conju- 
gate pairing rule of Lyapunov spectra for thermostated 
systems [3, pip, B, lo calized behaviors for Lyapunov vec- 
tors 0, H, 0, UU Il3l- an d the thermodynamic limit of 
Lyapunov spectra [lj, HiJ LL~J D3 an d so on have been 
discussed from this point of view. 

The stepwise structure of Lyapunov spectra is a typi- 
cal chaotic property of many-particle systems, which was 
found recently [16j- This stepwise structure appears in 
the Lyapunov exponents with smallest absolute value, 
and the dynamical structure of these Lyapunov expo- 
nents should reflect a slow and global behavior of the 
macroscopic system. Therefore, clarification of the step- 
wise structure of Lyapunov spectra (Lyapunov steps) is 
expected to make a bridge between the macroscopic sta- 
tistical theory and the microscopic chaotic dynamics. Re- 
markably, the Lyapunov steps accompany wavelike struc- 



tures in the associated Lyapunov vectors, namely the 
Lyapunov modes, which offer a useful tool to understand 
the origin of the stepwise structure of the Lyapunov spec- 
trum [17], [lj| [ljj, |2(j ■ Originally, these structures were 
observed in many-hard-disk systems, but very recently 
numerical evidence for the Lyapunov modes was reported 
for many-particle system with soft-core particle interac- 
tions [2l|. Some theoretical arguments have been pro- 
posed to explain this phenomenon, for exa mple, using 
random matrix theory |22tl23|. kinetic theory |24ll25ll26| . 
and periodic orbit theory |15|. etc. 

The first step to understand the Lyapunov step and 
mode is in the zero Lyapunov exponents and the as- 
sociated Lyapunov vectors. For instance, the two- 
dimensional system consisting of N particles with pe- 
riodic boundary conditions has six zero-Lyapunov expo- 
nents, and the associated Lyapunov vectors are linear 
combinations of the six basis vectors iV _1 / 2 (l, 0, 0, 0), 
iV-V^O, 1,0,0), N~ 1/2 (0, 0,1,0), N-^ 2 (0, 0,0,1), 
\p\~ 1 (Px,P v ,0,0), and |p| _1 (°, °,Px, P y ) where is a 
TV-dimensional null vector, 1 is an iV-dimensional vec- 
tor with all re-component equal to 1, and p = (p x ,Py) 
is the momentum vector with its ir-component p x and 
y-component p y . Here, the first (second) vector is 
associated with the translational invariance in the IE- 
direction (y-direction), the third (fourth) vector with 
the conservation of the x-component (y-component) of 
the total momentum, the fifth vector with the time- 
translational invariance (deterministic nature of the or- 
bit), and the last vector with the energy conservation 
[lj. This means that the six sets of the Lyapunov vec- 
tor components {Sx^jj, {Sy^jj, {Sp^jj, {Sp y f}j, 

{Sxf' /pxjjyp /p V j}j, and {Sp)$ IPxj,Wyi /Pyj}j can 
have equal components independent of the particle num- 
ber index j in zero-Lyapunov exponents. Here, we use 
the notation 8T^ = (<5q ( "\ <5p( n) ) for the Lyapunov 



vector corresponding to Lyapunov exponents X^ n ', and 
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• • • , Sp"^) are the spatial part and the momentum part 
of the Lyapunov vector, respectively. 

As the next step we regard the degeneracy of the zero- 
Lyapunov exponents and the structure of the correspond- 
ing Lyapunov vectors as the zero-th Lyapunov step and 
mode. This scenario was proposed first in Ref. |l8( , and 
was also discussed very recently in Ref. |20f. First of 
all, the Lyapunov steps for the two-dimensional system 
with periodic boundary conditions consist of two-point 
steps and four-point steps [171 Il8| , namely the number 
of Lyapunov exponents for one set of Lyapunov steps is 
six, which is equal to the number of the zero-Lyapunov 
exponents. It is also known that the stepwise structure 
of Lyapunov spectra is changed by violating the spatial 
translational invariance and the total momentum conser- 
vation, which also change the number of zero-Lyapunov 
exponents 18]. As a second point, some mode structures 
were observed in some of the above Lyapunov vector com- 
ponents, which should be constant in zero-Lyapunov ex- 
ponents. For example, a mode structure in the Lyapunov 
vector component Sy)j (the transverse spatial transla- 
tional invariance Lyapunov mode) is well known [17] . 
This mode is stationary in time, and appears in one of the 
two types of the Lyapunov steps. Ref. [lj| showed an- 
other mode structure in 5y\- /p y j, (the transverse time 
translational Lyapunov mode). This mode depends on 
time, and appears in other types of the Lyapunov steps. 
These Lyapunov modes are enough to categorize all the 
Lyapunov steps. Ref. [13) E3 a ^ so claim a moving mode 
structure in Sx n . 

However, there has not been enough evidence yet to 
confirm the above scenario for the Lyapunov steps and 
modes. For example, the mode structure in the momen- 
tum part of Lyapunov vectors has not been reported ex- 
plicitly. Besides, the phase relations of different modes, 
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for example, the modes in Sy^ /p y j and Sx, , have not 



been discussed. Another important problem is the time 
scale specified by a time-dependent Lyapunov mode, like 
the time-oscillation for the mode in Sy™ /p y j. The time- 
oscillating period is usually much longer than the mean 
free time of the system, and it should correspond to a 
collective movement, but quantitative evidence for it has 
not been shown explicitly. 

As an indicator for collective movements of many- 
particle systems, we can use the momentum auto- 
correlation functions, where collective movement may ap- 
pear as a time-oscillation behavi or [2711 as observed in 
many macroscopic models J23, |2|j, tM LiH |32| . The auto- 
correlation functions are accessible experimentally using 
neutron and light scattering techniques 32, 33, 34, 35]. 
As an essential aspect of auto-correlation functions is 
their role as response functions for the system. For ex- 
ample, linear response theory connects the time integral 



of the auto-correlation function with a transport coef- 
ficient 36]. However it should be emphasized that the 
auto-correlation functions themselves provide much more 
detailed information about the system. Linear response 
theory requires the time-integral and the thermodynamic 
limit of the auto-correlation function to calculate trans- 
port coefficients, and in this process information about 
short time-scales and a finite size effects are lost. For in- 
stance, the time-oscillation of the auto-correlation func- 
tions is one of the finite size effects, which linear response 
theory does not treat. Information on short time-scales 
and finite size effects in the auto-correlation functions 
also plays an important role in the generalized hydrody- 
namics [33 and generalized Fokker-Planck equation [38| . 

The purposes of this paper are two-fold. First we cal- 
culate all components of the Lyapunov vectors associated 
with the stepwise structure of the Lyapunov spectrum. 
They include the longitudinal and transverse components 
of both the spatial and momentum parts of the Lyapunov 
vectors. We demonstrate the wavelike structures in the 
components of the Lyapunov vectors, and specify their 
phase relations. These results support the above expla- 
nation for the origin of the Lyapunov steps and modes 
based on the zero-Lyapunov exponents and the associ- 
ated Lyapunov vectors. Spatial node structures of these 
Lyapunov modes are explained in terms of boundary con- 
ditions. It is emphasized that some of the Lyapunov 
modes show time-oscillating behaviors, and a particle 
number dependence in their time-oscillating periods. The 
second purpose of this paper is to discuss the connection 
between the time-oscillation of the Lyapunov modes and 
the momentum auto-correlation functions. Our central 
result is that the longest time-oscillating period of the 
Lyapunov modes is twice as long as that for the auto- 
correlation function for the longitudinal component of 
the momentum. This relation is satisfied irrespective of 
the number of particles in the system and the bound- 
ary conditions. We give a simple explanation for this 
relation. This result means that the time-oscillating be- 
havior of the Lyapunov vectors is reflected by a collec- 
tive movement within many-particle systems. This also 
gives some evidence to connect the Lyapunov mode, a 
tangent space property specified by Lyapunov vectors, 
to the auto-correlation function, which is a phase space 
property and is accessible experimentally. 

We use a quasi-one-dimensional many-hard-disk sys- 
tem, as this model allows fast numerical calculation of 
the Lyapunov exponents and vectors, and shows clear 
Lyapunov steps and modes. In general the numerical cal- 
culation of the Lyapunov spectrum and vectors are very 
time-consuming, and it is often difficult to get clear Lya- 
punov mode structures, particularly for time-dependent 
Lyapunov modes. If the above picture of the Lyapunov 
steps and modes based on universal properties such as 
the translational invariances and the conservation laws 
can be justified, then a simple model should be suffi- 
cient to convince us of their origin. This system also 
exhibits a clear oscillatory behavior in the longitudinal 



momentum auto-correlation function. A useful technique 
to get clear Lyapunov steps and modes is to use hard- 
wall boundary conditions. Although hard-wall bound- 
ary conditions destroy the spatial translational invari- 
ance and the total momentum conservation, and lead to 
different structure in the Lyapunov spectrum compared 
to periodic boundary conditions, it has been shown that 
there is a simple relation between the observed Lyapunov 
steps and modes and different boundary conditions jlSJ . 
Specifically, we use a quasi-one-dimensional system with 
hard-wall boundary conditions in the longitudinal direc- 
tion and periodic boundary conditions in the transverse 
direction. Usually, the hard-wall boundary conditions 
make numerical calculation slower compared to periodic 
boundary conditions, but in our system only the two par- 
ticles at each end of the system collide with the hard- walls 
so the effect is small. 

The outline of this paper is as follows. In Sec. [HJ the 
quasi-one-dimensional system is introduced. In Sec. IIIII 
we discuss the Lyapunov steps and modes. In Sec. HVI the 
momentum auto-correlation functions and their relation 
with the Lyapunov modes is discussed . Finally we give 
some conclusion and remarks in Sec. [V] 



II. QUASI-ONE-DIMENSIONAL SYSTEM 

The model considered in this paper is a quasi- 
one-dimensional many-hard-disk system. It is a two- 
dimensional rectangular system consisting of many hard- 
disks with the width of the system so narrow that the 
disk positions are not exchanged, thus the disks can be 
numbered from left to right. We assume that the mass 
m and the radius R of each disk is the same. In this case 
the quasi-one-dimcnsional system has a width L y that 
satisfies the condition 2R < L y < AR. Figure [I] gives a 
schematic illustration of the quasi-onc-dimensional sys- 
tem with the particles numbered 1, 2, • • • ,N (N: particle 
number) from the left to right. 



FIG. 1: A schematic illustration of a quasi-one-dimensional 
system used in this paper. The system shape is so narrow 
that particles always remain in the same order. Here, L x 
(L y ) is the length (width) of the system in the longitudinal 
(transverse) direction, and R is the radius of a particle. The 
dashed lines represent periodic boundary conditions and the 
solid lines represent hard-wall boundary conditions. The par- 
ticles are numbered 1, 2, • • ■ ,N from the left to the right. 



Originally, the quasi-one-dimcnsional many-hard-disk 
system was introduced as a system to easily and clearly 
observe the stepwise structure of Lyapunov spectrum 
(which we call the "Lyapunov steps") and the corre- 
sponding wave-like structure of Lyapunov vectors (the so 
called " Lyapunov modes" ) [lj| . Numerical observation of 



the Lyapunov steps and modes is very time-consuming, 
and even at present the Lyapunov spectra is limited to 
about 1000 particles [lTfl. Therefore it is valuable to ex- 
plore fast and effective ways to calculate them numeri- 
cally. It is well known that the rectangular system has a 
wider stepwise region than a square system with the same 
area and number of particles. This quasi-one-dimensional 
system is the most rectangular two-dimensional system 
possible. Ref. [lj| demonstrated that in the quasi-one- 
dimensional system we can clearly observe the structure 
of the Lyapunov modes. It is known that there are two 
kinds of Lyapunov modes: stationary modes and time- 
dependent modes. The stational Lyapunov modes are 
most easily observed because of their stable structure, 
but generally the observation of the time-dependent Lya- 
punov modes is much harder because of large fluctuations 
in their structure and their intrinsic time-dependence. 
The quasi-one-dimensional system is the first system 
in which time-oscillating Lyapunov modes were demon- 
strated [18J |. 

Another advantage of the quasi-one-dimensional sys- 
tem is that the particle interactions in this system are 
restricted to nearest-neighbor particles only, so we need 
much less effort to find colliding particle pairs numeri- 
cally compared with the fully two-dimensional system in 
which each particle can collide with any other particle. 
This leads to faster numerical calculation of the system 
dynamics. Besides, in this system, particle movement 
in the narrow direction is suppressed, compared to the 
longitudinal direction, and roughly speaking the particle 
sequence corresponds to the particle position. This leads 
to a much simpler representation of Lyapunov modes, 
which must be investigated as functions of spatial co- 
ordinates and time. One may also notice that in the 
quasi-one-dimensional system the system size is propor- 
tional to the particle number N, while for the square 
system it is proportional to y/~N. This implies that in 
the quasi-one-dimensional system, many-particle effects 
more evident than for the fully two-dimensional systems 
with the same number of particles. 

Another important point in the quasi-one-dimensional 
system is the effect of boundary conditions. Different 
from the square system, in which the boundary length 
is proportional to the square root of the system size, in 
the quasi-one-dimensional system the boundary length is 
proportional to the system size itself, therefore we can- 
not neglect its effect even in the thermodynamic limit. 
Actually Ref. 18] showed that the Lyapunov steps and 
modes depend strongly on boundary conditions. Bound- 
ary conditions change, not only the structure of Lya- 
punov steps and modes, but also the clearness of Lya- 
punov mode structure. For example, a system with 
purely hard-wall boundary conditions has no stationary 
Lyapunov mode and its corresponding Lyapunov steps 
and show much clearer time-oscillating Lyapunov modes, 
compared to a system with the periodic boundary con- 
ditions. Roughly speaking, hard-wall boundary condi- 
tions pin the positions of the nodes and thus lead to 



clearer Lyapunov modes. On the other hand, the numeri- 
cal calculation with the hard-wall boundary conditions is 
more time-consuming. This disadvantage is significant in 
the quasi-one-dimensional system with purely hard-wall 
boundary conditions. As another disadvantage of the 
system with purely hard- wall boundary conditions, we 
cannot investigate the stationary Lyapunov mode due to 
spatial translational invariance. As an optimal system 
we mostly consider a quasi-one-dimensional system with 
hard-wall boundary conditions in the longitudinal direc- 
tion and periodic boundary conditions in the transverse 
direction. We use the notation (H,P) for this boundary 
condition throughout this paper. In Fig. Qwe represent 
the boundary condition (H,P) as different types of lines 
on the boundaries: the bold solid lines signify hard-wall 
boundary conditions and the broken lines signify periodic 
boundary conditions. In this system we can get much 
clearer Lyapunov modes than for purely periodic bound- 
ary conditions and we observe both types of Lyapunov 
steps and modes. 

Although the quasi-one-dimensional many-hard-disk 
system with the boundary condition (H,P) may be ar- 
tificially introduced to investigate Lyapunov steps and 
modes in a fast and effective way, it is essential to note 
that the results from this system can be used to pre- 
dict Lyapunov steps and modes in more general systems, 
such as a fully two-dimensional system with purely peri- 
odic boundary conditions. Details of the relation of the 
Lyapunov steps and modes in quasi-one-dimensional sys- 
tems with different boundary conditions were given in 
Ref. [jj|- F° r example, the step widths of the Lyapunov 
spectrum (the spatial and time periods of the correspond- 
ing Lyapunov modes) in the system with the boundary 
condition (H,P) are halves (twice) the ones in the sys- 
tem with the purely periodic boundary conditions (P,P). 
It is also known that the structure of Lyapunov steps for 
the quasi-one-dimensional system is the same as the fully 
two-dimensional system. 

As we discussed above, the main reason to use the 
quasi-one-dimensional system with the boundary condi- 
tion (H,P) is to get clear Lyapunov steps and modes in a 
fast numerical calculation. On the other hand, these ad- 
vantages may not assist the calculation of the momentum 
auto-correlation function. In the quasi-one-dimensional 
system, collisions of a particle are restricted to its two 
nearest-neighbor particles only, so it may be supposed 
that specific types of collisions like the "back scattering 
effect" play an important role in this system. The back 
scattering effect, which comes from a reversal of the ve- 
locity of a particle by a collision with the nearest neighbor 
particle, can lead to a negative region of the momentum 
auto-correlation function. One may also suppose that 
in the quasi-one-dimensional system a collective motion 
may be enhanced, because the movement of particles in 
the narrow direction is very restricted. This may lead to a 
clear time-oscillation of the momentum auto-correlation 
function, as will be actually observed in Sect. IIVI of 
this paper. It may also be noted that boundary condi- 



tion effects on the momentum auto-correlation function 
are not well known. Because we need to know about 
them to be able to guess the relation between the auto- 
correlation functions and Lyapunov modes in different 
boundary conditions, we will discuss boundary condition 
effects on momentum auto-correlation function briefly in 
Sect. ITV71 

In this paper we use units where the mass m and the 
particle radius R are 1, and the total energy E is N 
(except in Sec. IIIIDfl . For the numerical calculations, 
the system lengths are chosen as L x — 1.5NL y + 2R and 
L y = 2R(1 + 1CP 6 ) for the quasi-one-dimensional system 
with the boundary condition (H,P). 



III. LYAPUNOV STEPS AND MODES 

In this section we discuss the Lyapunov steps and 
modes in the quasi-one-dimensional system with bound- 
ary condition (H,P). Part of these results have already 
been presented in Ref. jl8|. and here we complete this 
presentation. Some of the discussions omitted in Ref. 
[18f were the relation between the time-oscillating Lya- 
punov mode proportional to the momentum and the lon- 
gitudinal Lyapunov modes, the Lyapunov modes for the 
momentum parts of Lyapunov vectors, and the particle 
number dependence of the time-oscillating period of the 
Lyapunov modes. 

We introduce the Lyapunov spectrum as the ordered 
set of the Lyapunov exponents X^ n \ n — 1,2,- • ■ , 4/V, 
where A^ > \& > • • ■ > X^ N) . The notations Sq^ = 

(SXj,8yj) and Sp" = (Sp^ , 6p£ ) are used for the 
spatial and momentum components, respectively, of the 
j'-th particle Lyapunov vector corresponding to the n-th 
Lyapunov exponent A^™' . 

For a numerical calculation of the Lyapunov spectrum 
and the Lyapunov vectors, we used the numerical algo- 
rithm developed by Benettin et al. | 39| and Shimada 
et al. ;40j (Also See Refs. HH E2). This algorithm is 
characterized by regular re-orthogonalizations and renor- 
malization of set of Lyapunov vectors, which can be done 
after each particle collision in a many-hard-disk system. 
Usually the Lyapunov steps and modes appear after a 
long trajectory calculation, and we typically calculated 
trajectories of more than 5 x 10 5 particle collisions to get 
the Lyapunov spectra and vectors shown here. 

The main purpose in this section is to investigate the 
time-oscillating structures in the transverse and longitu- 
dinal Lyapunov modes. As we mentioned in Sect. ITTlthc 
numerical calculation of Lyapunov spectra and Lyapunov 
vectors for many-particle systems is time-consuming even 
in the quasi-one-dimensional system, and in this pa- 
per we present results of 100 particle systems. In such 
a rather small system the structure of the Lyapunov 
modes has large fluctuations in space and time, which 
prevents the appearance of clear mode structures. An- 
other problem appears in the investigation of the Lya- 



punov modes proportional to the momentum. We in- 
vestigate these modes by calculating such the quantities 

5x n /p x j, which should show a time-oscillating wave- like 
structure, if the Lyapunov vector components has modes 
proportional to momentum. However, fluctuations of this 
quantity may be enhanced in cases where the momen- 
tum component has a small absolute value. To overcome 
these technical problems, and to visualize this structure 
as clearly as possible, we take a local time-average of the 
Lyapunov modes. More concretely, for investigations of 
the modes of the Lyapunov vector components Sx? and 

Sy\ we take their arithmetic average over 8N collisions 
using data just after collisions, and plot them as func- 
tions of the same local time arithmetic average of Xj and 
the first collision number of time interval taking the lo- 
cal time average. For modes in the quantities Sx ' " /p X j 

fyj /Pyji SPyj I Pxj, and dp]™ /pyj we take the same lo- 
cal time average, except that if the absolute value \p X j \ 
(\p y j\) of the momentum component is less than 10% of 
the averaged momentum amplitude y / 2ME/N then we 
exclude the data at that time from the local time average. 
We use the notation (•••)* for such a local time average 
in this paper. 



A. Lyapunov steps 

Figure [21 is the stepwise structure of the Lyapunov 
spectrum normalized by the largest Lyapunov exponent 
for the quasi-one-dimensional many-hard-disk system of 
100 particles with (H,P) boundary condition. The entire 
positive branch of the Lyapunov spectrum is shown in the 
inset to this figure. In a Hamiltonian system, the neg- 
ative branch of the Lyapunov spectrum takes the same 
absolute value as the positive branch of the Lyapunov 
spectrum from the conjugate pairing rule: \( iN - n + 1 ) — 
— \( n \ n = 1, 2, • • • , 2N [43| . so they are omitted in Fig. 

121 



FIG. 2: Stepwise structure of the Lyapunov spectrum nor- 
malized by the largest Lyapunov exponent for the quasi-one- 
dimensional system with hard-wall boundary conditions in 
x-direction and periodic boundary conditions in y-direction. 
The particle number is N = 100. The white-filled (black- 
filled) circles correspond to the stational (time-oscillating) 
Lyapunov modes. Inset: The full spectrum of the positive 
branch of the normalized Lyapunov spectrum. 



wall boundary condition in the x-direction. Half of these 
4 zero-Lyapunov exponents appear in Fig. [21 

The stepwise structure of the Lyapunov spectrum in 
this system is one and two-point steps. These two kinds 
of Lyapunov steps accompany different mode structures 
in the Lyapunov vectors: one is the stationary modes, 

and the other is time- 



as discussed in subsections IIII Bl 
oscillating modes, as discussed in the subsection IIII CI 
Here we count the sequence of Lyapunov steps from the 
zero Lyapunov exponents, so A^ 200 ^ and A*- 199 -* are the 
zero exponents, A^ 198 ' is the first one-point step, A 1 - 197 ) 
and A^ 196 ) are the first two-point step, A^ 195 -* is the second 
one-point step, and A^ 194 ' and A^ 193 - 1 are the second two- 
point step, see Fig. [21 



B. Stationary Lyapunov modes 

First we discuss the Lyapunov mode corresponding to 
the first and second one-point steps in Fig. Fig. |3| 
shows the graph of the Lyapunov vector components cor- 
responding to the first and second one-point step as a 
function of the collision number n t and the normalized 
local time average (xj}t/L x of the x-component of the 
particle position. Here, Fig. |3{a) is for the local time- 



averages of (5; 



„(198)\ 



and (Sy 



(198) v 



for the first one-point 



step A^ 198 ), and Fig. E|L>) is a similar graph for the sec- 
ond one-point step A^ 195 ^. These one-point steps are in- 
dicated by arrows in Fig. [2] Both graphs have the same 
collision number interval [524000,569600]. On the base 
of each of Figs. EUa) and (b) we give contour-plots of the 



,(198) \ 



,(195)\ 



transversal modes {Sy, )t and {Syy a '''}t, respectively, 
in which the dotted lines, the solid lines, and the bro- 
ken lines correspond to the levels -0.08, 0, and +0.08, 
respectively. 



FIG. 3: Local time averages of (Sxj)t and {5yj)t for the 
first two-point step (exponent 198) in (a), and for the second 
two-point step (exponent 195) in (b). These are shown as 
functions of the collision number n t and the normalized local 
time average {xj)t/L x of the x-component of the j-th parti- 
cles position (with the system length L x ) . The corresponding 
Lyapunov exponents A^ 198 ' and A^ 195 ' are indicated by arrows 
in Fig. |21 The base of each graph are the contour plots of 
the transverse Lyapunov modes at the levels —0.08 (broken 
fines), (solid lines), and +0.08 (dotted lines). 



This system has 4 zero-Lyapunov exponents, which 
come from the conservation of the y-component of the 
total momentum and the center of mass, energy conser- 
vation, and the deterministic nature of the orbit. Note 
that the ^-component of the total momentum and the 
center of mass are not conserved, because of the hard- 



In Fig. El we recognize spatial wave-like structures in 
the transverse components of {Sy? )t f° r n — 198 and 
195 which are stationary in time (at least in a time inter- 
val of more than 45 x 10 3 collisions as shown in this fig- 
ure) . These wave- like structures are very nicely fitted by 
sinusoidal functions .18]. Note that in the calculation of 
Lyapunov vectors we used the numerical algorithm which 



imposed renormalization of the Lyapunov vectors at ev- 
ery collision, so that the amplitudes of any component of 
the Lyapunov vector must be less than 1. It should be 
emphasized that anti-nodes in the modes appear at the 
end of the system in the x-direction. By comparison the 



»\ 



for 



amplitudes of the longitudinal components of (8: 

n = 198 and 195 are extremely small. These observations 

suggest that the Lyapunov mode corresponding to the k- 

th one-point Lyapunov step is approximately represented 

by 
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M*0) 



m;, cos 1 —X.j 



(1) 

(2) 



j = 1,2, ••• ,N, corresponding to the Lyapunov expo- 
nents A"^*", in which a/b is a constant and fj,(k) is the 
Lyapunov index corresponding to the fc-th one-point step 
of the Lyapunov spectrum. Here we take the origin of the 
x-component of the spatial coordinate to be Xj — 0, so 
that an ambiguity in spatial phase can be removed in Eq. 

©■ 



C. Time-oscillating Lyapunov modes 

Now we discuss the remaining Lyapunov modes, which 
correspond to the two-point steps of the Lyapunov spec- 
trum. Figure 01 shows the graphs of the local time 



„(197) 



(197) 



,(197)\ 



averages of {8x) ')t, {Sx) ' /p x j)t, {Sy) )t, and 
(fy) /Pyj)t as functions of the collision number n* and 
normalized local time averaged position (xj)t/L x of the 
j-th particle. The number of particles is N = 100, and 
the four graphs are for the same collision number interval 
n t e [535200,569600]. These correspond to the first ex- 
ponent on the first two-point step, exponent 197, in Fig. 
121 In these figures, we draw a contour plot on the base 
of each three dimensional graph at levels —0.08 (dotted 
lines), (solid lines), and +0.08 (broken lines). 



FIG. 4: Local time averages of (8xj)t, {8xj/p x j)t, (8yj)t, and 
{Syj/p y j)t for the Lyapunov modes corresponding to the first 
exponent of the first two-point step (A' 197 '), as functions of 
the collision number nt and the normalized local time aver- 
age {xj)t/L x of the i-component of the position of the j-th 
particle. The base of each graph is a contour plot of the three 
dimensional graph at the levels —0.08 (dotted lines), (solid 
lines), and +0.08 (broken lines). 



We can easily recognize spatial wave-like structures 
with time-oscillations in Figs. Ufa) and (d). In Fig. Ufa), 
the longitudinal Lyapunov vector component (8x) )t 
has nodes at the ends of the quasi-one-dimensional sys- 
tem, and the wave-length is given by 2L X . On the other 



hand, in Fig. EJd), the transverse Lyapunov vector com- 
ponent (Sy, /Pyj)t has anti-nodes at the end of the 
system and has a node at the middle, although its wave- 
length is given by 2L X as for the longitudinal mode of 
Fig. Ha). _ 

There is also a time-oscillating wave-like struc- 
ture in the longitudinal Lyapunov vector component 

(5xj ' /p X j)t, as shown in Fig. Efb). This structure 
is different from the one shown in Fig. Ufa) associated 
with the same longitudinal Lyapunov vector component 
Sx ■ . It has anti-nodes at the ends of the system, and 
has a node in the middle. Its wave-length is 2L X . These 
characteristics suggest that although there are large fluc- 
tuations in the middle of the system, the time-oscillating 
wave-like structure in Fig. 0fb) is the same as that in Fig. 
BJcl). In Fig. 0Jc) it is rather difficult to recognize any 
structure. Roughly speaking, it is just random fluctua- 
tions, in middle of the system but the amplitude of such 
fluctuations is small compared to region at the end of the 
system. However such small amplitude fluctuations are 
required for consistency with Fig. EJd), namely the fact 
that in this region the value of Sy) I Pyj is small, so 

the value of Sy :■ itself should be small with an almost 
position independent momentum p y j . 

Next we discuss the phase relations for the Lyapunov 
modes of the first two-point steps, corresponding to ex- 
ponents A'- 196 -' and A 1 - 197 ) (the black-filled circles with 
brace underneath in Fig. [2J . Figure shows the con- 
tour plots of the local time averages of (8xj)t, (Sxj/p X j) t , 
(Syj/p y j)t, and for A^ 196 ) and A^ 197 ^ as functions of the 
collision number n t and the normalized local time aver- 
age (xj) t /L x of the position of the j-th particle at levels 
—0.08, 0, and +0.08, for a 100-particle system. The six 
graphs in Fig. [SJhave the same collision number interval 
[535200,569600], and Figs. EJa), (c), and (e) correspond 
to Figs. Qla), (b), and (d), respectively. 



FIG. 5: Contour plots of the local time averages of {8xj)t, 
(8xj/p x j)t, (8yj/Pyj)t, for the first two-point steps (A*- 197 ' 
(a,c,e) and A' 196 ' (b,d,f)), as functions of the collision num- 
ber nt and the normalized local time average (xj)t/L x of the 
position of the j-th particle in the same collision number in- 
terval [535200, 569600]. The dotted lines, the solid lines, and 
the broken lines are contour lines at the levels —0.08, 0, and 
+0.08, respectively. 



Figure shows that the two Lyapunov exponents for 
the same two-point step have the same structure of Lya- 
punov modes, but they are orthogonal in time, namely 
node lines of the Lyapunov modes corresponding to the 
Lyapunov exponent A^ 197 -* correspond to anti-node lines 
of the ones of A^ 196 ^ . We also notice that the node lines of 
the Lyapunov modes in (Sx)" /p x j)t and (Sy, \ /Pyj)t co- 
incide with each other in space and time (n = 197, 196), 



on the other hand the Lyapunov modes in (Sx n )t and 

(Sx " /p x j)t are orthogonal in space and time at the same 
Lyapunov index. 

The above discussions based on Figs. 0] and [3] (and 
similar observations of Lyapunov modes in the other two- 



point steps of the Lyapunov spectrum) lead to the con- 
jecture that the spatial part of Lyapunov vector compo- 
nents Sxj and Sy" corresponding to the Lya- 
punov exponents constructing the /c-th two-point step are 
approximately expressed as 
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1, 2, • • • , N with constants a' k , a 



k' "fei u fe) "fc 



and /%. It ber of particles N. 



should be noted that large fluctuations in the Lyapunov 
mode represented in middle of Figs0Ib),[3c), andHJd) 
can come from the second terms on the right-hand sides 
of Eqs. 10 and Q}. On the other hand the effect of the 
first terms on the right-hand side of Eqs. J2J) and (@J 
does not appear explicitly in FigsEJa), OJa) andCJb), 
because the factor p X j in these terms distributes their 
contributions randomly and these terms disappear after 
taking local time averages. 



FIG. 



Amplitudes 
~ 3) )t (circles), 



Oti I and 

(2JV-3) 



«1 
/Pxj)t 



of Lyapunov 
(triangles), and 



modes (Sx- )t (circles), (5x 

(Sy\ ~ /Pyj)t (squares) as functions of ^j2mE/N in the 



(2JV-3) 

quasi-one-dimensional system of 50 hard-disks with (H,P) 
boundary condition. The broken line is a fit of the amplitude 
|ai| to a constant function y — £ with a fitting parameter £, 
and the solid line is given by y = (,x. 



D. Energy dependence of Lyapunov mode 
amplitudes 

In the expressions ©, ®, © and © for the time- 
oscillating Lyapunov modes, the quantities a' k , a' k , a' k 
and a' k are introduced simply as coefficients of the lin- 
ear combination of the longitudinal spatial translational 
invariance Lyapunov mode and the time translational 
invariance Lyapunov mode. However it is important 
to note that these coefficients are related to each other 
through the normalization of the Lyapunov mode. 

We consider the condition that the Lyapunov 



mode vector (Sx^ 
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v' = v(k),u(k) — 1 is normalizable. 
leads to the approximate relations 
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(7) 



y/2mE/N 

141 - 141, (8) 

141 ~ 141 (9) 

with the mass m(= 1), the total energy E and the num- 



In Fig. El we show the amplitudes \a[\ and |5i|, which 
are obtained by fitting the Lyapunov modes (Sx 



(S. 



„(2Af-3) 



/p x j)t, and (Sy 



(2N-3) 



(2N- 
3 



Iv 



'vj/t 



to sinusoidal func- 



tions multiplied by constants, as functions of y / 2mE/N. 
To get the data for this figure we used the quasi-one- 
dimensional system of 50 hard-disks with (H,P) bound- 
ary conditions, and calculated the amplitudes \a'i\ and 
| fill for different total energies E. In Fig. we fitted 



the amplitude |fij| for the mode (S: 



,(2AT-3)\ 



(circles) to 



a constant function y — £ (the broken line) with a fitting 
parameter value £ « 0.179, and the solid line is given by 



This y — S, x using this value of £. The amplitudes |a^| for the 



modes (& 



„(2JV-3) 



IVxj)t (triangles), and (Sy 



(2JV-3) 



Iv 



'volt 



(squares) are reasonably on the line y = £x, and these 
results support the relation JJJ, and also suggest that the 



,(2JV-3)v 



is independent 

.(2JV-3) 



amplitude |fi' fc | for the mode (8: 

of ^/2mE/N. The amplitude |a' fe | for (Sxf N ~' S) /p xj ) 

(triangles) and (Sy^ ~ /Pyj)t (squares) in Fig. al- 
most coincide with each other, and it gives support to 
the claim that the coefficients a' k (in the first term) on 



the right-hand side of Eqs. JHJ and JHJ coincide. 

It may be noted that the normalization procedure in 
the Benettin's algorithm, which we used to calculate the 
Lyapunov exponents and vectors in this paper, requires 
the normalization of Lyapunov vectors including both 
their spatial part and momentum part, but in the above 
argument we assumed the normalizability of the spatial 
part only of the Lyapunov vectors. This can be justified 
by the fact that as shown in Sec. IIII Fl the spatial part 
and momentum part of Lyapunov vectors constructing 
Lyapunov modes show almost the same mode structure, 
so each of them should be independently normalizable. 



results explain the spatial node structures shown in Fig. 
[3] The spatial node structure of the stationary Lyapunov 
modes in Syj corresponding to the one-point steps can be 
explained in this way. Because Syj varies sinusoidally and 
must satisfy the reflection symmetry, it must be either a 
node (if the sign changes) or an anti-node (if the sign is 
invariant). Hence, in this case the Lyapunov mode in Syj 
should have an anti-node at the hard-walls. 



Lyapunov modes in momentum components of 
Lyapunov vectors 



E. Spatial node structures of the Lyapunov modes 
and reflections in the hard-walls 



The spatial node structure of the Lyapunov modes can 
be explained using the collision rule for particles with 
hard- walls. 

For (H,P) boundary condition the particle collisions 
with the hard-walls in the ^-direction cause a change 
in the sign of the x-component of the momentum with 
the remaining components of the phase space vector un- 
changed: 
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(10) 

(11) 

(12) 
(13) 

Similarly, in this type of collision the x-components of the 
Lyapunov vector change their signs while the remaining 
components unchanged: 



(14) 
(15) 
(16) 
(17) 



Note that in the x components of Lyapunov vector Sxj 
changes its sign as well as 5p X j, which is different from 
the phase space vector. 

The important point is that a system with hard-wall 
boundaries is equivalent to an infinite system generated 
by reflecting the positions and velocities of all particles 
(in the hard wall) and by changing the signs of all x- 
components of the Lyapunov vectors at the hard-wall. 
That is explicitly incorporating the reflection symmetries 
for the phase space vector and the Lyapunov vector at 
hard walls. If the modes of the entire system are con- 
nected smoothly sinusoidal functions at the hard walls, 
then this condition requires that the mode for the quan- 
tity Sxj has a node at a hard wall, because it changes 
sign there. On the other hand, the quantities Sxj/p X j 
and Syj / Pyj do not change their signs at hard walls, so 
these modes should have anti- nodes at hard walls. These 



So far, we have discussed only the spatial components 
of the Lyapunov vectors. In this subsection we discuss 
briefly the Lyapunov modes appearing in the momentum 
parts of Lyapunov vectors. 

One of the few differences between spatial and momen- 
tum components of Lyapunov vectors is that the am- 
plitudes of the momentum components are often much 
smaller than those of the corresponding spatial compo- 
nents [ll| . This makes Lyapunov modes for the momen- 
tum parts of the Lyapunov vectors less clear than the cor- 
responding spatial components. However, basically the 
structure of the Lyapunov mode for the momentum part 
of the Lyapunov vector is quite similar to the correspond- 
ing spatial component. For this reason, in this subsec- 
tion we omit a detailed discussion of the phase relations 
of multiple Lyapunov modes for the momentum parts of 
Lyapunov vectors, and just show that there are certain 
modes structures in the momentum components of Lya- 
punov vectors corresponding to the Lyapunov steps. 



FIG. 7: Global time averages (Sp™.) for the Lyapunov expo- 
nents corresponding to the first, second and third one-point 
step of the Lyapunov spectrum [n = 198 (circles), 195 (trian- 
gles) and 192 (squares)] as functions of the normalized global 
time average {xj)/L x of the position of the j-th. 



Figure shows mode structure of Sp y ™ corresponding 
to the first three one-point steps (n = 198,195 and 192) as 
functions of the normalized particle positions. This struc- 
ture is stationary in time, so we took their global time 
average over 200iV collisions, using the notation (■ ■ • ) for 
this global time-average without the suffix t. These mode 
structures are stationary in time, and are similar to the 
ones for the corresponding spatial part Sy^ discussed in 
Sect. iHLBl 

Figure [S] shows contour plots of time-oscillating 



J197) 



Lyapunov modes for (Sp xj ') t , (Sp xj /p x j)t, and 
(Sp , /p y j)t as functions of the collision number rit and 
the normalized local time average (xj)t/L x in the first 
two-point step. We used the same collision number in- 
terval [535200, 569600] in Fig. Has in Fig. The mode 



FIG. 8: Contour plots of the local time averages of (5p x j)t, 
{Spxj/Pxj)t, and {8pyj/p y j)t for the first exponent (197) of the 
first two-point step as a function of the collision number nt 
and the normalized position of the j-th particle (in the same 
collision number interval [535200, 569600]). The dotted lines, 
solid lines, and broken lines are contour lines at the levels 
-0.0018, 0, and +0.0018, respectively. 



FIG. 9: Particle number (TV) dependence of the period Ti ya of 
the time-oscillations of (Sx- ~ )t (circles), (5x, ~ /p x j)t 

(triangles), and (8y- ~ /p y j)t (squares) in the quasi-one- 
dimensional system with (H,P) boundary condition. The data 
is fitted to the function y = a + j3x 2 with the fitting parame- 
ters a and /3. The inset: Particle number dependence of the 
mean free time r with a fitting function y = 7/2 with the 
fitting parameter 7. 



structures in Figs. El (a), (b) and (c) are almost the same 
as Figs. EI a), (c) and (e) for the corresponding spatial 

components (Sxy 97 ') tl (5xf 97 > /p XJ ) t , and (Syf 97 ' /p y j)t, 
respectively, although their oscillating amplitudes are 
much smaller than those of the corresponding spatial 
components. 

The spatial mode structures of the momentum com- 
ponents of Lyapunov vectors are explained by the same 
reflection property at hard- walls, which was discussed in 
the previous subsection 1111 El 



G. Particle number dependence of the oscillating 
periods 



shown in Fig. gfa), (b), and (d) for N = 100. Fig. 
El shows that the three time-oscillations associated with 



(Sx} '}t, (Sx} '/p X j)u and (Sy} /p y j)t all have the 
same period. In Fig. [5|the data is fitted to a quadratic 
function y — a + (3x 2 with the fitting parameter values 
a ss 17.9 and /3rj 1.65 The inset to this figure shows the 
mean free time r as a function of the number of particles 
N. The N dependence of r is nicely fitted to the function 
y = j/x with the value 7 sw 1.91 of the fitting parameter 
7. Noting that the period, in real time, is given by rT\ ya 
approximately, these results suggest that the period of 
the Lyapunov modes is proportional to the number of 
particles TV, namely the system size. 



In Sect. IIII CI we showed that the quantities (Sx 
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(Sx j /p X j)t, and (5yj /p y j)t corresponding to the two- 
point Lyapunov steps show time-oscillating behavior. 
Now we consider how the time-oscillating period of those 
Lyapunov modes depends on the number of particles N 
for the quasi-one-dimensional system. 

We evaluate the collision number interval for the time- 
oscillation of Lyapunov modes as follows. As shown in 
the proceeding subsection IIII CI the Lyapunov modes re- 
lated to the quantity (8x™ )t (the quantities (5x " /p X j)t 

and (Syj 1 /p y j)t) have an anti-node in the middle (at the 
end) of the system in the x-direction (n = 2N — 3 and 
2N — 4) for (H,P) boundary condition. Using this prop- 
erties we took 6 data points for the quantity (Sx " }t (the 

quantities (Sxj /p x j)t and (Syy' /p y j)t) («• = 2N—2) in 
the middle (at the end) of the system with (H,P) bound- 
ary condition. These data are fitted to a sinusoidal func- 
tion y — asa\{(2'Kx/Tiy a ) + b} with fitting parameters 
a, b, and Ti yai which leads to a numerical estimation of 
the period TJ yo of the time-oscillation of the Lyapunov 
modes. The collision number interval Ti ya can be trans- 
lated into a real time interval by multiplying by the mean 
free time r, if necessary. 

Figure is the graph of the period Ti ya of the 

time-oscillations of (5x)j )t (circles), (Sx^ /p X j)t (tri- 

,(197), 



3 
angles), and (Sy { - /p y j)t (squares) in the quasi-one- 
dimensional system with (H,P) boundary condition, as 
functions of the number of particles N. Spatial and 
temporal behavior of these quantities has already been 



FIG. 10: The quantity L x / (TT ac f) as a function of the num- 
ber of particles N. Here, Ti ya is the collision number interval 
given by the average of the three collision number intervals 
Tiya for time-oscillations of the Lyapunov vector components 
(5xf N ~ 3) )t, (5xf N ~ 3) /p X j)t, and (5yf h ' /p y j)t in the first 
two point step of the Lyapunov spectra. The line is given by 
y — 1, which is the thermal velocity yj E/(MN). 



Now, we investigate the time-oscillating period of the 
Lyapunov modes in a different way. Figure 1101 shows 
the quantity L x /(rT ac f) using the system length L Xl 
the mean free time r and the averaged time-oscillating 
periods rT ac f of the Lyapunov modes in (Sx, )t, 



(2N-3) 



(2N-3) 



(Sxj /p X j)t, and (8yf" °' /p y j)t in the first two- 

point step of the Lyapunov spectrum. This figure sug- 
gests that this quantity is almost 1 independent of the 
particle number N, therefore equal to the thermal veloc- 
ity y/E/(MN). 



IV. AUTO-CORRELATION FUNCTIONS 

In this section we discuss another property of the quasi- 
one-dimensional system, namely the time-oscillation be- 
havior of the momentum auto-correlation function. This 
is a typical measure of the collective behavior of many- 
particle systems. We connect this behavior with the 



10 



time-oscillating behavior of the Lyapvmov modes, sug- 
gesting that the time-oscillation of the Lyapunov modes 
can also be regarded as a collective mode. 

We calculate numerically the auto-correlation func- 
tions C v (t) for the 77-component of the momentum based 
on the normalized expression C v (t) = (7^(t)/C' 7 ,(0), in 
which C v (t) is defined by 



C v (t) 



lim 



1 



t^+oo (N 2 -N t + 1)T 

N 2 r T 
X 



N 2 ^T 

^2 / ds Pvj( s + t )Pvj( s )^ ( 18 ) 
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77 = X or y. Eq. (|18|l includes a time- average and an 
average over some of particles (from the Ni-th particle 
to the 7V2-th particle) in the middle of the system. (Note 
that number the particles 1, 2, • • • ,N from left to right 
in the system, as shown in Fig. Q]) In actual calculations 
we choose A^i = [(JV + l)/2]-5 and N 2 = [(JV+l)/2] + 5 
with [x] as the integer part of the real number x. This 
means that we take into account only 11 particles in 
the middle of the system in the calculation of the auto- 
correlation function C v (t). (In this paper we consider 
the case N > 40 > 11.) It should be noted that us- 
ing (H,P) boundary condition the auto-correlation func- 
tion for particles near hard-walls are different from the 
ones for particles in the middle of the system, as dis- 
cussed in Appendix^ Especially, the momentum auto- 
correlation function of particles near hard-walls do not 
show clear time-oscillating behavior. To get the clearest 
time-oscillating behavior for the auto-correlation func- 
tion C v (t) and to get less hard- wall boundary condition 
effects, we exclude the auto-correlation functions of par- 
ticles near hard- walls in the calculation of C v (t) . 

If the system is ergodic, the value of the auto- 
correlation function (|18fl will be independent of the ini- 
tial condition. To get the results for the auto-correlation 
function in this section we take a time-average of the 
auto-correlation function over more than 2 x 10 6 colli- 
sions. For convenience, in figures the auto-correlation 
functions are shown as functions of the correlation num- 
ber n t , but it can always be translated into the real time 
t by multiplying by the mean free time r. 



of the momentum auto-correlation functions. In this sys- 
tem the mean free time is given by t w 0.0188. From 
Fig. 1111 it is clear that the momentum auto-correlation 
function has a strong direction dependence and shows a 
time-oscillating behavior in C x . 



FIG. 11: The auto-correlation function C x and C y for the x- 
and y-components of momentum, respectively, as functions 
of the collision number nt- Main figure: Linear- linear plots 
of C x and C y as functions of nt- The inset: Log-log plots 
of the absolute values \C X \ and \C y \ of the auto-correlation 
functions as functions of tit- Here, the broken line is a fit of 
the graph of \C X \ to an exponential function, and the line is 
a fit of the graph of \C y \ to a ^-exponential function (defined 
byEq. J2"Hh 



In the beginning, the auto-correlation function C x for 
the x-component of momentum decays exponentially in 
time. To show this point, in the inset to Fig. 1111 we 
fitted the beginning of the graph of \C X \ to an exponential 
function 



Gi(z) = exp{-a:'a;} (19) 

with the fitting parameter value a' sa 0.0385. 



FIG. 12: The time-oscillating part of auto-correlation func- 
tion C x for the a;-component of momentum as functions of the 
collision number n t . Here, the line is the fit to the product of 
a sinusoidal and an exponential function. 



The significant point about the auto-correlation func- 
tion C x is its time-oscillating behavior. To show this 
behavior explicitly we show Fig. 1121 as an enlarged graph 
of the time-oscillating part of C x , which is already shown 
in Fig. ^] This time-oscillation accompanies a time de- 
cay, so we fitted this graph to the product of a sinusoidal 
and an exponential function G2(x), namely 



A. Momentum auto-correlation functions and their 
direction dependence 

Figure 1111 contains the momentum auto-correlation 
functions C x and C y for the momentum components in 
the x- and y-directions, respectively, as a function of 
the collision number nt in the quasi-one-dimcnsional sys- 
tem with N = 100 and (H,P) boundary condition. The 
main figure in Fig. 1111 is a linear-linear plot of the auto- 
correlation functions C x and C y , while its inset is a log- 
log plot of the graph of the absolute values \C X \ and \C y \ 



G 2 {x) = A e~ p x sin 
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T, 



ucj 



-x + t 



(20) 



with fitting parameters A, f3' , T ac f, and £. The time 
oscillating part of the auto correlation function C' x is 
nicely fitted to this function with the parameter values 
A « 0.0209, /?' w 5.17 x 10~ 5 , T ac} w 8.29 x 10 3 , and 
£ ~ 1.62. This also gives us a way of numerically eval- 
uating the oscillation period T ac f of the auto-correlation 
function C x . We note the quasi-one-dimensional system 
shows a much clearer time-oscillating behavior of the mo- 
mentum auto-correlation function than a fully two- (or 
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three-) dimensional system. One may ask whether the 
damping behavior of the envelope of time-oscillation of 
C x is best fitted to a power-law function, like the slow 
damping of the long time behavior of C y , rather than to 
an exponential function as assumed in Eq. (|20l) . (Actu- 
ally the data in Fig. ^]is not sufficient to decide between 
exponential and power decay.) This point is discussed 
further in Appendix IbI 

On the other hand, the auto-correlation function C v 
for the y-componcnt of momentum shows a significantly 
different behavior from C x . This comes from the specific 
shape of the system and the boundary condition. As 
shown in Fig. 1111 the damping of the momentum auto- 
correlation function C y in the y-direction is much slower 
than in the cc-direction. This is simply explained by the 
fact that in the quasi-one-dimensional system particle 
collisions occur mostly as head-on collisions, so that the 
change in the y-component of momentum in a collision 
can be much smaller than the change of the x-component 
of momentum. In the inset to Fig. 1111 the graph of 
the auto-correlation function for the y-component of mo- 
mentum is fitted to the " re-exponential function" F K (x), 
which is defined by 



F K (x)= y/l + (a" rex) 2 



a kx 



1/k 



(21) 



in x > with fitting parameters a" and re. In the colli- 
sion number region shown in Fig. Illl the auto-correlation 
function C y in the y-direction is positive, so this fitting 
can be for C y as well as for \C y \. From the definition, in 
the limit as re = the function F K (x) becomes the ex- 
ponential function: lim K ^o F K (x) — exp{— a"x}, noting 
F„(0) = 1 and dF K (x)/dx = -a"F K (x)/y/l + (a" K x) 2 , 
so this function is a one-parameter deformation of the 
exponential function [44J . The important properties of 
this function are that it is approximated by an exponen- 
tial function at small x and is approximately a power 
function in a large x. 



F K (x) 



(2na"x)~ 1/K 



in x « 1 

in x » 1 and na" x > 



(22) 



The fitting of the numerical data for the auto-correlation 
function C y to the re-exponential function with parame- 
ter values a" « 0.00358 and re w 1.44 is very satisfactory, 
and this implies that this auto correlation function decays 
exponentially in the beginning, (like the auto-correlation 
function in the x-direction), and decays as a power func- 
tion after that, at least in the time scale shown in Fig. ITT1 
(This does not mean that the auto-correlation function 
C y decays as a re-exponential function in any time-scale. 
See Appendix [B] about the auto-correlation function C y 
at much longer time-scales than shown in Fig. UTT ) 



B. Particle number dependence of the 

auto-correlation function and a relation with the 

time-oscillation of the Lyapunov modes 

We have shown the two kinds of time-oscillation be- 
haviors in the quasi-one-dimensional system: one for the 
Lyapunov mode and another for the momentum auto- 
correlation function. Now we show numerical evidence 
to connect these two behaviors. 



FIG. 13: The period Ti ya (in collision numbers) of the time- 
oscillation of the Lyapunov mode as a function of the pe- 
riod of time-oscillation of the longitudinal momentum auto- 
correlation function T ac f. Data points are obtained from nu- 
merical calculation of the quasi-one-dimensional system with 
(H,P) boundary condition for different numbers of particles 
N - 40, 50, 60, • • • , 100. The line is given by the function 

J lya ^-L acf • 



Fig. E| is the graph of the time-oscillating period 
Ti ya of the Lyapunov modes as a function of the time- 
oscillating period T ac f of the momentum auto-correlation 
function C x . Here, the time-oscillating periods are calcu- 
lated for N = 40, 50, 60, • • • , 100, and the time-oscillating 
periods Ti ya are calculated as the average of the colli- 
sion number interval Ti ya for time-oscillations of the Lya- 



.(2JV-3)v 



(<5 



(2JV-3) 



111 



'xj/U 



punov vector components (5; 

and (Sy. /Pyj)t in the first two-point step of the Lya- 

punov spectra. (As shown in Sect. IIII Gl these three 
oscillating periods Ti ya of the Lyapunov modes take al- 
most the same values.) In Fig. ^|we also show the line 
given by the function Ti ya = 2T ac f. The numerical data 
for the time-oscillations in Fig. ^|is nicely fitted to the 
function, 



-L l y a ~ ^ acf • 

This is the main result of this paper. 



(23) 



In Table [I] we summarize, not only the values of the 
time-oscillating periods Ti ya r and T ac fT of the Lyapunov 
modes and the momentum auto-correlation function in 
real time, but also the data about the iV-dependences of 
the damping properties of the auto-corrclation functions 
C x and C y . They include the mean free time r, the 
exponential damping times r/a' and t//3' (for C x ) and 
r/a" (for C y ), and the power 1/re of the damping of C y at 
long time. Here, values of a', (3' and T ac f (a" and re) are 
derived by fitting the auto-correlation function C x (C y ) 
to Eqs. (O and (gOJl (Eq. (|2TJ>)- From this table it is 
clear that the exponential damping times r/a' and r/a" , 
and the power 1/re are almost independent of the particle 
number N. On the other hand, the exponential damping 
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TABLE I: Time-oscillating periods and decay factors for the Lyapunov modes and the momentum auto-correlation functions. 
Here, N is the number of particles, r is the mean free time, and Ti ya is the collision number interval given by taking an 
average of the three collision number intervals Ti ya of time-oscillations of the Lyapunov vector components (8xj)t, (&Xj/p Xj )t, 
and {Syj/p y j)t corresponding to the first two point step of the Lyapunov spectra. The parameter a' is given by fitting the 
beginning of the longitudinal momentum auto-correlation function C x as a function of the collision number nt to an exponential 
function y = exp{— ax}. The parameters and T ac f are given by fitting the time-oscillating part of the same function C x to 
the function y = Aexp{—/3'x} sin[(27r/T ac /):z;-|-£]. The parameters a" and k are given by fitting the transverse momentum auto- 
correlation function C y as a function of the collision number nt to the K-exponential function y — [^/l + (a"nx) 2 — a" nx] 1 ' K 
(Eq. ED ). 





Lyapunov 


mode 




Cx 




Cy 


N 


r 


Tly a T 


rjd 


r//3' 


T ac f T 


r/a" 


1/k 


40 


0.0489 


123.5 


0.493 


74.2 


63.0 


5.33 


0.732 


50 


0.0380 


154.7 


0.483 


102.5 


76.0 


5.26 


0.716 


60 


0.0326 


197.3 


0.505 


156.1 


96.6 


5.37 


0.675 


70 


0.0275 


223.1 


0.493 


196.1 


107.9 


5.26 


0.608 


80 


0.0238 


256.7 


0.487 


204.5 


125.0 


4.74 


0.577 


90 


0.0210 


283.3 


0.487 


252.0 


139.9 


4.89 


0.624 


100 


0.0188 


305.7 


0.488 


362.7 


155.5 


5.27 


0.693 



time t/(3' of the time oscillation of the auto-correlation 
C x increases as N increases. (We have already discussed 
the iV-dependence of the mean free time r in Sec. IIIIGI ) 



C. Boundary condition effects 

So far, we have concentrated into the quasi-one- 
dimensional system with hard-wall boundary conditions 
in the ^-direction and periodic boundary conditions in 
the y-direction, namely the (H,P) boundary condition, 
for technical convenience in the anal ysis of the Lyapunov 
modes. On the other hand, in Ref. [lg we have already 
discussed and compared Lyapunov steps and modes in 
the different boundary conditions: the purely periodic 
boundary conditions (P,P), the purely hard- wall bound- 
ary conditions (H,H), and periodic boundary conditions 
in the x-direction and hard-wall boundary conditions in 
the y-direction (P,H) as well as the boundary condi- 
tion (H,P). In this section we carry out a similar discus- 
sion for the momentum auto-correlation functions C x in 
these different boundary conditions. Figure. 1141 contains 
schematic illustrations of these boundary conditions. 

For meaningful comparisons between the different 
boundary conditions we use the same mass m and ra- 
dius R for the particles, and the same number of parti- 
cles (N = 50). Using the set of the lengths (L x ,L y ) to 
define the size of the system in the x and y directions for 
(H,P) boundary conditions, then we use (L x — 2R, L y ) 
for (P,P) boundary conditions (L x — 2R, L y + 2R) for 
(P,H) boundary conditions, and (L x , L y + 2R) for (H,H) 
boundary conditions. This gives the same effective area 
for particles to move in each of the four systems. This 
also means that the mean free time r in these four types 
of the boundary conditions will be the same (Concrete 
numerical values of r are given in Table [H] ) . 



FIG. 14: Schematic illustrations of the four boundary con- 
ditions (P,P), (P,H), (H,P) and (H,H) used in quasi-one- 
dimensional systems. Here, (P,P) is the purely periodic 
boundary conditions, (P,H) is periodic boundary conditions 
in the x-direction and hard-wall boundary conditions in the 
{/-direction, (H,P) is hard-wall boundary conditions in the in- 
direction and periodic boundary conditions in the y-direction, 
and (H,H) is the purely hard-wall boundary conditions. The 
dashed lines and the solid lines on the boundaries represent 
periodic boundary conditions and hard-wall boundary condi- 
tions, respectively. 



FIG. 15: Auto-correlation functions C x for the ^-component 
of the momenta as functions of the collision number n t for the 
boundary conditions (P,P), (P,H), (H,P) and (H,H). The sys- 
tems are quasi-one-dimensional systems consisting of 50 hard- 
disks. We observe: (a) Exponential decay region in the begin- 
ning of the damping of the auto-correlation function C x as a 
linear-log plot. The dotted line and the broken line are the fits 
for the cases (P,P) and (H,P) and the cases (P,H) and (H,H) 
to exponential functions, respectively, (b) Time-oscillating 
region of the auto-correlation functions C x as a linear-linear 
plot. The four graphs of the auto-correlation functions C x are 
fitted to the functions y — A exp{— fH'x} sin{(27r /T ac f)x + £} 
with the fitting parameters A, /3', T ac f and £. 



Fig. 1 151 shows the auto-correlation functions C x for the 
x-component of the momenta in quasi-one-dimensional 
systems consisting of 50 hard-disks with boundary con- 
ditions (P,P), (P,H), (H,P) and (H,H) as functions of the 
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collision number fit- Here, Fig. I15f a~) is the beginning 
part of the autocorrelation functions C x , and is given 
as a linear-log plot to show their exponential decay as 
straight lines. In this figure the fits to the exponential 
function (|19|l with the fitting parameter a' are given 
for the cases (P,P) and (H,P) and the cases (P,H) and 
(H,H) separately. The dotted line is the fit for the 
cases (P,P) and (H,P) with the fitting parameter values 
a' rs 0.0765, and the broken line is for the cases (P,H) 
and (H,H) with the fitting parameter values a' « 0.0597. 
Fig- H5f b^l is the time-oscillating part of C x in the 
four different boundary conditions. In this figure, each 
auto-correlation function is fitted to a function 1)20(1 with 
the fitting parameters A, 13' , T ac f, and £. The values 
of these fitting parameters are {A, /3',T ac /,^) « 
(0.0422,0.00108,1.00 x 10 3 , 1.37) for " (P,P), 
(A,{3',T acf ,£) sw (0.0447,0.000803,1.24 x 10 3 ,1.33) 
for (P,H), (A,(3',T ac f,£) w (0.0398,0.000369,2.04 x 
10 3 ,1.56) for (H,P), and (A,(3\T acf ,0 » 

(0.0403, 0.000246, 2.53 x 10 3 , 1.56) for (H,H). 

Figure lTBT a - ) shows that the boundary condition in the 
y-direction has a strong effect on the auto-correlation 
function C x even at short time. In the cases (P,P) and 
(H,P) the auto-correlation function C x shows faster expo- 
nential decay than for (P,H) and (H,H), as shown in the 
difference of the value of the fitting parameter a' for the 
exponential fitting function (|19fl . In the collision number 
region shown in Fig. H5f aL the effects of the bound- 
ary conditions in the x-direction in the auto-correlation 
function C x appear after showing their initial exponen- 
tial decays, and the autocorrelation function C x for (P,P) 
(and (P,H)) decays faster than C x for (H,P) (and (H,H)). 

On the other hand, Fig. I15f b~) shows that in all the 
boundary conditions the auto-correlation functions C x 
show time-oscillations, but with different oscillating pe- 
riods T ac f. In this figure we can recognize that the second 
peak of C x for (P,P) coincides with the first peak of C x 
for (H,P). A similar coincidence appears in the second 
peak of C x for (P,H) and the first peak of C x in (H,H). 
Actually the fitting parameter value of T ac f for (P,P) 
((P.H)) is approximately the half the value of T ac f for 
(H,P) ((H,H)). This can be simply explained by the fact 
that replacing the periodic boundary conditions with the 
hard-wall boundary conditions, it needs twice the time 
for a particle perturbation to come back to the same po- 
sition. 



TABLE II: The mean free time r, the averages of the longest 
time-oscillating time period Ti ya r for the Lyapunov vector, 
and the time-oscillating time period T ac fT of the longitudi- 
nal momentum auto-correlation function C x for the differ- 
ent boundary conditions (P,P), (P,H), (H,P) and (H,H) in 
a quasi-one-dimensional system consisting of 50 hard-disks. 
Tiy a is the average over three collision number intervals of 
Ti ya . Here, Ti ya is the average of the time-oscillations of the 

three components (5:rA fc) } t , {dXj /p X j)t, and {Sy^ /p y j)t for 
the first Lyapunov step which has time-oscillating Lyapunov 
modes (ft = 2iV-5 for (P,P), ft = 2N -2 for (P,H), ft = 27V-3 
for (H,P), and k = 2N - 1 for (H,H)). 



Boundary 


T 


Tl ya T 


T acf r 


(PP) 


0.0369 


37.4 


77.0 


(P,H) 


0.0371 


45.8 


91.4 


(H,P) 


0.0380 


77.3 


154.5 


(H,H) 


0.0383 


96.8 


194.5 



for the quantities (6: 



.(*)\ 



(5x™/ P: 



xj)t, and {5yf ] /p 



yj/t- 



The Lyapunov indices k are chosen from the Lyapunov 
exponents in the first Lyapunov step which has time- 
oscillating behavior of its Lyapunov modes. Our result 
supports the conjecture that the relation (|23() is satis- 
fied for all boundary conditions (P.P), (P,H), (H,P) and 
(H,H). 



D. An explanation for the relation of 

time-oscillation periods of the Lyapunov mode and 

the momentum auto-correlation function 

As we have shown, the relation Ti ya = 2T ac f (Eq. (|23|) ~) 
between the largest time-oscillating period Ti ya of the 
Lyapunov modes and the time-oscillating period T ac f of 
the momentum auto-correlation function is independent 
of the number of particles N and the boundary condi- 
tions. In this subsection we discuss a possible explana- 
tion for this relation, which is a physical argument rather 
than a strict mathematical proof. 

We consider a momentum component p x (t), like 
the x-component of the momentum in the quasi-one- 
dimensional system, which shows a time-oscillating be- 
havior in its auto-correlation function with a frequency 
w ac /: 



Finally we show a relation between the time-oscillating 
periods of the Lyapunov mode and the momentum auto- 
correlation function for different boundary conditions. 
In Table \n\ we summarize the mean free time r, the 
time-oscillating time period Ti ya r of the largest Lya- 
punov mode, and the the time-oscillating time period 
Tacf''' of the momentum auto-correlation function C x for 
the four kinds of the boundary conditions (P,P), (P,H), 
(H,P) and (H,H). Here, the period fi ya is evaluated as the 
arithmetic average of the collision number intervals Ti ya 



p x (t)* p x (0) ~ <j>{t) e 1 ' 



(24) 



where we use the notation X(t)*X(0) for the auto- 
correlation function for any complex quantity X(t) with 
complex conjugate X(t)* . Here <j)(t) is the damping enve- 
lope of the auto-correlation function p x (t)* p x (0), which 
can be an exponential decay, <fi(t) ~ exp{— at} with a 
positive constant a (See Appendix iBlh 

On the other hand, we represent the momentum pro- 
portional and time-oscillating term in the Lyapunov vec- 
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tor, like the first terms on the right-hand sides of Eqs. 
© and 10, as 



Sq x ~ il)i(t)p x (t)e l 



(25) 



with a frequency u>i ya , where tp\(i) is the decay enve- 
lope of the amplitude of 5q x > an( i m &y show an expo- 
nential divergence (or contraction) following the corre- 
sponding Lyapunov exponent. Now, we assume that if 
the quantity dq x oscillates persistently in time, then its 
auto-correlation function 5q x {t)* 8q x (0) should oscillate 
in time with the same frequency LOi yai namely 



Sq x (t)* Sq x (0) ~ ^(tje*"-* 

with a new envelope function "02 (i)- 

It follows from Eqs. JSJ, J23) and (J2SJ that 



(26) 



Mt)e 



i^lyat 



Vi (ty mo) Px(ty p x (p) e-™""* 
-- M*Y Mo) W) e* UaB '~ u,va)t , 

which immediately leads to 



Mt) - Mt)*M0)cb(t) 

Ulya ~ W ac //2. 



(27) 
(28) 



The time-oscillating periods T ac f and r; ya of the momen- 
tum auto-correlation function and the Lyapunov mode 
are given by T acf — 27r/(rw ac/ ) and T tya ~ 2n/ (rwj„ ). 
Using this point and Eq. I|28|l we obtain our relation (I23|) . 
Note that the above explanation for Ti ya = 2T ac f is inde- 
pendent of the number of particles N and the boundary 
conditions. 



(2iV-3) 



FIG. 16: The normalized auto-correlation function C t 

for the longitudinal Lyapunov vector component 8xa as 

a function of the collision number rit . The system is a quasi- 
one-dimensional system consisting 50 particles with boundary 
condition (H,P). The numerical data is well fitted to a sinu- 
soidal function multiplied by an exponential function. 



In the above explanation, the assumption i|26[) is cru- 
cial, so it may be useful to demonstrate this behavior 
numerically. Figure 1161 shows that the auto-correlation 

(IN 3) 

function C,„„ „ for the longitudinal Lyapunov vec- 



y lya, X 

tor component Sx 



(2N-3) 



normalized by its initial value 



(about 0.0203) in an quasi-one-dimensional system of 
50 particles with (H,P) boundary condition. In the 

auto-correlation function C\ ~ its mean value is sub- 
tracted, and an average over the auto-correlation func- 
tions of 11 particles in the middle of the system is taken. 



Here, the Lyapunov index n = 2N — 3 of the Lya- 
punov vector component 5x™ is chosen so that the cor- 
responding Lyapunov step is the first two-point step as- 
sociated with a time-oscillating Lyapunov mode. In Fig. 
1161 the numerical data is fitted to a sinusoidal function 
multiplied by an exponential function, namely the func- 
tion Q20JI. with the fitting parameter values A sw 0.967, 
/?' « 1.03 x 10~ 5 , T acf = Tacf « 4.12 x 10 3 , and £ « 1.54. 
This time-oscillating period T ac f for the auto-correlation 
function for the longitudinal Lyapunov vector component 
coincides almost exactly with the time-oscillating period 
Tacf ~ 4.07 x 10 3 of the corresponding Lyapunov mode. 
This coincidence of the time-oscillating periods supports 
our assumption i|26*|) 4_5(. 



V. CONCLUSION AND REMARKS 

In this paper we have discussed the relation between 
the wave-like structure of Lyapunov vectors and the time- 
oscillating behavior of the momentum auto-correlation 
functions in quasi-one-dimensional many-hard-disk sys- 
tems. The quasi-one-dimensional system is a narrow rect- 
angular system in which the x components of the particle 
positions remained in the same order. This system was 
proposed as a many-particle system which shows clear 
stepwise structure of the Lyapunov spectrum (the Lya- 
punov steps) and wave-like structure of the associated 
Lyapunov vectors (the Lyapunov modes). Using this sys- 
tem, we showed that there are two types of Lyapunov 
modes in the spatial and momentum components of the 
Lyapunov vectors corresponding to the two kinds of steps 
in the Lyapunov spectrum: one is stationary in time and 
the other involves a time-oscillation. Here, the time- 
oscillating Lyapunov vectors consist of a simple time- 
oscillating part plus a momentum proportional time- 
oscillating part in the longitudinal components, while 
the transverse time-oscillating Lyapunov vectors consist 
of a momentum proportional time-oscillating part only. 
We revealed the phase relation for these time-oscillating 
Lyapunov modes. It was shown that the system length 
divided by the time-oscillating period of the Lyapunov 
modes is independent of the number of particles at the 
same density, and is of order of the thermal velocity. 
After discussing these wave-like structures of the Lya- 
punov vectors, we connected them to the time-oscillation 
of the momentum auto-correlation. The time-oscillation 
of the auto-correlation function appears in the longitu- 
dinal component of the momentum, and its envelope de- 
cays exponentially in time. The main point is that the 
largest time-oscillating period of the first time-oscillating 
Lyapunov modes is twice as long as the time-oscillating 
period of the momentum auto-correlation function. We 
showed that this relation is independent of the number of 
particles and the boundary conditions (constructed from 
combinations of periodic and hard-wall boundary condi- 
tions) . A simple explanation is given for this relation. It 
was also shown that the auto-correlation function for the 
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transverse component of the momentum is nicely fitted 
to the K-exponential function, implying that it decays ex- 
ponentially at the beginning and decays as a power after 
that. 

In this paper we considered mainly a specific boundary 
condition for the quasi-one-dimcnsional system: (H,P) 
hard-wall boundary conditions in the longitudinal di- 
rection and periodic boundary conditions in the trans- 
verse direction. The system with this boundary condi- 
tion exhibit a much clear wavelike-like structure of Lya- 
punov modes than the purely periodic boundary condi- 
tions (P,P), which is a big advantage for quantitative dis- 
cussions of the Lyapunov modes. Using the (H,P) bound- 
ary condition, the spatial translational invariance in the 
longitudinal direction is violated, and it leads to a differ- 
ent Lyapunov step structure and auto-correlation func- 
tions, compared with the (P,P) boundary conditions. For 
example, in (H,P) the step widths of the Lyapunov spec- 
trum are half of the ones in (P,P), and individual particles 
can have different momentum auto-correlation functions 
due to the back scattering effect of the hard-wall (See 
Appendix II V (Jj l while the momentum auto-correlation 
function is particle-independent for the (P,P) boundary 
condition. However, as discussed in Rcf. 18] for the Lya- 
punov modes and in Sec. IIVCI for the auto-correlation 
functions, there is a simple relation connecting the re- 
sults obtained from different boundary conditions, so we 
can predict some of the results of the other boundary 
conditions from the results for (H,P). 

The mode structure of Lyapunov vectors discussed in 
this paper is related to the structure of the Lyapunov 
vectors associated with zero-Lyapunov exponents. As ex- 
plained in the introduction of this paper, there are sets 
of Lyapunov vector components which take a constant 
value independent of the particle index, and that these 
quantities corresponding to the stepwise structure of the 
Lyapunov spectrum have wavelike structures. These are 
connected with the spatial and time translational invari- 
ances and the energy and momentum conservation laws. 
However we need to be careful when making a connec- 
tion between the conservation laws (or the translation 
invariances) and the Lyapunov modes. For example, in 
a system with hard-wall boundary conditions the spa- 
tial translational invariance is violated, but even in such 
systems the mode structure in the Lyapunov vector com- 
ponent Sxj (or 6y\ ■ ), can be observed. However, a 
scenario which suggests that translational invariance is 
only evident when it is observed in the zero-Lyapunov 
exponent modes, will not predict these observed longitu- 
dinal modes. 

It should be noted that a time-dependence of the Lya- 
punov modes may not always appear as a time-oscillating 
behavior. Ref. [ljj claim that the spatial wave of the 
Lyapunov vector " moves" at a specific speed in the square 
system consisting of many hard-disks. It is interesting 
to know how these different behaviors, one oscillating in 
time and another moving with a speed, can appear. 

In some papers, an understanding of the Lyapunov 



modes was attempted based on an analogy with the 
hydrodynamic modes |19|. Actually, in both cases the 
conservation laws like the total momentum conservation 
and the energy conservation play an essential role, and 
the longitudinal mode shows a time-dependent behav- 
ior. However it is important to know that the deter- 
ministic nature of orbits also plays one of the essential 
roles in the Lyapunov modes and leads to momentum- 
proportional time-oscillating components of Lyapunov 
vectors, although such a characteristic does not appear 
explicitly in the hydrodynamic mode. In this sense, it is 
still an open question to see how hydrodynamic modes, 
which have no concept of a phase space trajectory, can 
incorporate time translational invariance. 

From results of this paper, it is suggested that there is a 
connection between existence of the stepwise structure of 
Lyapunov spectra and the time-oscillations of momentum 
auto-correlation functions. It is well known that the step- 
wise structure of the Lyapunov spectra appears clearly in 
rectangular systems rather than in square systems at the 
same density. It is possible to get a similar result for 
the time-oscillation of the momentum auto-correlation 
function? For example, in a square system with a small 
number of hard-disks we cannot observe the stepwise 
structure of the Lyapunov spectrum, and in this case the 
time-oscillation of the momentum auto-correlation func- 
tion does not appear. Therefore, the time-oscillations 
of the auto-correlation function may be useful to under- 
stand the condition for existence of the Lyapunov steps 
and modes. In this sense, for example, it may be in- 
teresting to investigate the time-correlation function in 
systems with soft-core particle interactions in which the 
observation of the Lyapunov steps is much harder, and 
less direct than in systems with hard-core interactions. 
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APPENDIX A: MOMENTUM 

AUTO-CORRELATION FUNCTION OF 

INDIVIDUAL PARTICLE 

In this appendix we discuss the momentum auto- 
correlation function of individual particles in the quasi- 
one-dimensional system with hard-wall boundary condi- 
tions in the longitudinal direction and periodic bound- 
ary conditions in the transverse direction (the boundary 
condition (H,P)). Different from purely periodic bound- 
ary conditions, the hard-wall boundary in (H,P) violates 
the translation invariance in the longitudinal x-direction, 
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and this implies that different particles can have different 
momentum auto-correlation functions. 



We introduce the auto-correlation function 



JJ) 



(*) of 



FIG. 18: Exponential decay region of the auto-correlation 



the j-th particle in the ^-direction at time t (J — 
1, 2, • • • , N and 7/ = x, y) based on the normalized ex- 
pression c^\t) = c^\t)/c^\0) with c\p(t) defined by 



„( 2 ) 



_ functions c x (circles), c x (triangles), and c x (squares) 



(r,0) 



rU) 



C^(t) 



lim — 

T^+ooT 



for the x-components of momenta, as functions of the col- 
lision number nt on a linear-log plot. The solid line (the 
dotted-broken lines) is a fit of the graph of c x (the graphs of 
c x and c x ) to an exponential function. Similar graphs are 
given for the auto-correlation functions for the y-component 



(i) 



,(2) 



ds Pkj(s + t)pkj{s). (Al) of momenta: c y ' (inverted triangles), c y ' (pluses), and c 



„(50) 



Using this quantity ctf (t), the auto-correlation function 
C v (i) defined by Eq. l(T%|l is simply given by C n (t) — 



-A' 2 



7&), 



[i/(n 2 -n 1 + i)]j:^ Ni c v 



(t). In this appendix we 

show graphs of c^' as a function of the collision number 
n t pa t/r in the quasi-one-dimensional system consisting 
of 100 hard-disks. We number the particles 1, 2, • • • , N 
from the left to right, as shown in Fig. ^ so i f° r example, 
the 1st and N-th. particles are closest to the hard-walls. 



FIG. 17: Auto-correlation functions for c x , c x and c x 
for the ^-components of momenta of the 1st, 20-th, and 50-th 
particle, respectively, as functions of the collision number nt- 
The system is a quasi-one-dimensional system of 100 hard- 
disks with (H,P) boundary condition. The inset: Enlarged 
graphs in the small magnitude part of the auto-correlation 
functions. 



The first important point of the individual auto- 
correlation functions is that the time-oscillation in the 
x-direction is weak for particles near the walls. This is 
shown in Fig. El m which we plot the auto-correlation 
function c x for the j-th. particle in the x-direction for 
j = 1 (the nearest particle to the left hard- wall), j = 50 
(the particle most distant from the hard- walls), and 
j = 20 (a particle between these two extremes). In 
Fig. El the main figure is the full data for these auto- 
correlation functions, and its inset is an enlarged graph 
to emphasize the time-oscillating part. This figure shows 
that we cannot recognize a time-oscillating behavior in 
the auto-correlation function c x of the particle nearest 
to the hard-wall, although a clear time-oscillation can be 

recognized in c x of the particle in a middle of the sys- 
tem. We can see a time-oscillating behavior in the auto- 
correlation function c x , but its amplitude is smaller 
The positions of the nodes of the 



than that of c x 



time-oscillations of c x and c x almost coincide with 
each other. 

Another difference between individual particle auto- 
correlation functions appears in the short time scale. Fig- 
ure 1181 shows the momentum auto-correlation functions 

(i) ^~ 
eff , j — 1,2 and 50, as functions of the collision num- 
ber rit showing the initial damping behavior (77 = x,y). 



(crosses). The broken line (the dotted line) is a fit of the 



„C0 



(2) 



„( 5 °h 



graph of c y (the graphs of c y and c y ) to an exponential 
function. 



These auto-correlation functions show an exponential de- 
cay, which we present as a linear-log plot (straight lines 
imply exponential decay). This figure shows that the x- 
component (y-component) of the auto-correlation func- 
tion of the particle nearest the hard-wall decays faster 
(slower) than the ones of other particles, while in any of 
them the damping behavior is nicely fitted to an expo- 
nential function. (In Fig. El the graphs are fitted to the 
exponential function (|19|) with the fitting parameter val- 
ues a' k, 0.0457 for c x 1} (solid line), a' w 0.0369 for c x 2) 
and ci 50) (dotted-broken line), a' w 0.00184 for c { y ] (bro- 
ken line), and a' sw 0.00368 for c y ' and Cy ' (dotted line) 
). This difference may come from the different types of 
collisions experienced. For the particle nearest the wall, 
half of the collisions will be with the wall and the other 



half with the neighboring particle. The x-component of 
the momentum is drastically changed (namely it changes 
the sign of p x j)> so it may cause a faster decay of the au- 
tocorrelation function in the x-direction than for other 
particles. On the other hand, collisions with the wall 
effect the y-component of the momentum much less, be- 
cause it is invariant under wall collisions, and does not 
cause loss of memory of p y j , and this can also explain a 
slower decay of the momentum auto-correlation for the 
1st (and iV-th) particle in the y-direction compared to 
other particles. 



FIG. 19: The negative region of auto-correlation functions 
Cx , Cx , Cx , and c x for the ^-components of momenta, 
as functions of the collision number n t as a linear-linear plot. 
The inset: Absolute values of the same auto-correlation func- 
tions as functions of n t as a log- log plot. 



Another point of difference of the auto-correlation 
functions of individual particles is a negative region which 
appears after their initial exponential decay. It may be 
meaningful to mention that a negative region of mo- 
mentum auto-correlation function has drawn attention 



in some previous works J46J, |4jJ |4£| • To discuss such a 
negative region, in Fig. 1191 we show the collision number 

n t dependence of the auto-correlation functions c x , c x , 
Cx , and c x .To emphasize such a negative region of the 
auto-correlation function, we show a log-log plot of the 
absolute values of the same quantities as functions of rit 
in the inset to Fig. [2UJ As shown in this figure, a negative 
region of the auto-correlation functions appears after the 
initial exponential decay and before the time-oscillates 
appear. The collision number (or time) at the bottom 
of this negative region of the auto-correlation function 
increases, and the amplitude of the bottom decreases, 
as the particle is further from the hard-wall (namely as 

(1) 
the particle index j in c x increases from j — 1 to 10 

in Fig. I20|) . This phenomenon can be explained by the 
backscattering effect of the hard- walls. Such a back scat- 
tering effect is stronger (so the amplitude of the negative 
region is stronger) in a particle closer to a hard-wall, As 
well, the time interval to react to the presence of the wall 
is longer (so the time at the bottom of the negative re- 
gion is later) in a particle far from the hard-wall. This 
kind of behavior is not observed in a system in which the 
boundary conditions in the x-direction are periodic. 



FIG. 20: The region before the start of the time-oscillation 



of auto-correlation functions for c x ' . <■ 
functions of the collision number nt 



(10) „(20) „(35) 



and c. 



(.-.()) 
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APPENDIX B: DAMPING BEHAVIOR OF 

MOMENTUM AUTO-CORRELATION 
FUNCTION IN A LONG TIME INTERVAL 

In this appendix we discuss two points about the mo- 
mentum auto-correlation function in the long time inter- 
val: (i) The shape of the envelop of the time-oscillation 
in the auto-correlation function C x for the x-component 
of the momentum, and (ii) The behavior of the auto- 
correlation function C' y of the y-component of the mo- 
mentum on a much longer time scale than that shown in 
the text of this paper. 



FIG. 21: (a) Absolute values \C X \ and \C y \ of the 

auto-correlation functions for the x-component and the y- 
component of the momentum, respectively, as functions of 
the collision number n t presented as a log-log plot. The solid 
line is a fit of the envelope of the time-oscillating part of the 
auto-correlation function C x to an exponential function. The 
dotted line is a fit of the auto-correlation function C y to a 
K-exponential function (Eq. 12H - ) and the dotted broken line 
is a fit to an exponential function in the region where there 
is a deviation from the /^-exponential function, (b) The time- 
oscillating part of the auto-correlation part as a graph of the 
absolute value \G X \ as a function of the collision number nt 
presented as a linear- log plot. The broken line is a fit to a 
sinusoidal function multiplied by an exponential decay func- 
tion (Eq. 120|l i. and the solid line is its envelope, which is the 
same as the solid line in Fig. (a). The system is a quasi-one- 
dimensional system consisting of 50 hard-disks with (H,P) 
boundary condition. 



After such a negative region of the auto-correlation 
function the time-oscillating part appears. Figure |2"UI 
shows the collision number n t dependence of the auto- 
correlation functions for c x , c x , c4 and c x in the 
collision number region before the time-oscillation of the 
auto-correlation starts (about nt ~ 6000 in Fig. I20fl . 
The negative peak of the auto-correlation function (dis- 
cussed in the previous paragraph and indicated by the 
arrows in Fig. 12011 moves to a longer collision number rit 
as the particle index j increase from j = 10 to 35 in Fig. 
|2*U1 On the other hand, the time oscillation of the auto- 
correlation function starts from about n t ~ 6000 which 
is independent of the particle index, although the am- 
plitude of the time-oscillation is large for a particle far 
from the hard- walls. Moreover the time-oscillating pe- 
riod of the auto-correlation function is almost indepen- 
dent of the particle index. These characteristics of the 
time-oscillation of the auto-correlation function suggest 
that the time-oscillating behavior of the auto-correlation 
function reflects a collective movement of the system. 



Figure [2] shows the absolute values of x and y auto- 
correlation functions \C X \ (\C y \) as functions of the colli- 
sion number nt in a quasi-one-dimensional system of 50 
hard-disks with (H,P) boundary condition. In Fig. 12 If a) 
these graphs are plotted as log- log plots, while in Fig. 
I21f bi the graph for \C X \ is plotted as a linear-log plot. 
The collision number interval in this figure is about ten 
times as long as the previous ones in this paper, and we 
took a much longer time-average (eg. over 10 9 collisions) 
to get this data. 

As shown in Sec. IIV Al the auto-correlation function 
C x for the x-component of the momentum decays ex- 
ponentially initially. After the initial decay, the time- 
oscillating region of C x starts. We fitted this region of 
C x to a sinusoidal function multiplied by an exponen- 
tial function, namely Eq. (|20|l . with the fitting parame- 
ters A w 0.0402, 13' sw 0.000368, T ac} s=s 2.03 x 10 3 and 
£ rs 1.53 as the broken lines in Fig. I21f bi. The solid lines 
in Fig. EflT al and (b) are the envelope y = Aexp{—f3'x} 
of this function. In order to see its exponential behavior 
we show, in Fig. I21f bb the linear- log plot of \C X \ for the 
time-oscillating region of C x , in which the exponential 
decay is represented as a straight line. In this linear-log 
plot the local maximum points of \C X \ are clearly on a 
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straight line. 

In Sec. IIV Al we also showed that the auto-correlation 
function C y for the y-component of the momentum is 
nicely fitted to a K-exponential function (|21l) . This is also 
shown in Fig. I21f a) as the fit line to the K-exponential 
function with fitting parameter values a" ~ 0.00746 and 
k « 1.48 (dotted line). However, Fig. I21f a) shows that 
there is a deviation from this functional form on a long 
time scale. Such a deviation is significant in the colli- 



sion number region nt > 10000 in this graph. We fitted 
the auto-correlation function C y to an exponential func- 
tion y = A' exp{— a'"x} with fitting parameter values 
A' « 0.0369 and a"' « 6.25 x 10~ 5 (the dotted-broken 
in Fig. 12 If all in the region where C y deviates from the 
K-exponential. 
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